Enhanced Pinning of Vortices in Thin Film Superconductors by Magnetic Dot Arrays 
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We study the pinning of vortices in thin film superconductors by magnetic dots in the London 
approximation. A single dot is in general able to pin multiple field-induced vortices, up to a sat- 
uration number n s , which can be much larger than one. However, the magnetic field of the dot 
also creates intrinsic vortices and anti- vortices, which must be accounted for. In a ferromagnetic 
dot array, the intrinsic anti-vortices are pinned only interstitially. Much stronger pinning effect is 
expected of an antiferromagnetic dot array. Possible realizations of various magnetic configurations 
are discussed. 

PACS numbers: 74.60. Ge, 74.80. Dm, 75.50.Tt 

Artificially patterned sub-micron magnetic structures 
hold a great promise not only for magnetic device and 
storage technology, but also as a tool of fundamental re- 
search when used as a means of controlling other physi- 
cal systems, such as the two-dimensional electron gas Q 
and vortices in superconductors It is an experi- 

mental fact that when a type-II superconducting film is 
deposited over a regular array of magnetic dots, resis- 
tivity of the film exhibits a series of sharp minima as a 
function of the applied magnetic field, with the positions 
of the minima being integer multiples of the geometrical 
"matching field" of the magnetic lattice [0. In contrast, 
no periodic pinning was found in a similar system with 
a regular array of non-magnetic defects [^|. This sug- 
gests the importance of the magnetization of the dots 
in producing low resistivity, although the exact flux pin- 
ning mechanism has not been fully understood. A recent 
study by Lyuksyutov and Pokrovsky Jsj explored vari- 
ous statistical mechanics issues that might arise from the 
complexity of the magnet-superconductor interaction. In 
this work, we focus on the pinning mechanism itself, the 
understanding of which will lead us to desirable magnetic 
dot structures that enhance vortex pinning for a range of 
applied magnetic fields. 

We study the low temperature properties of a super- 
conducting thin film such as Nb deposited on top of a 
regular array of magnetic dots, separated by a thin in- 
sulating layer to suppress the proximity effect (Fig. 1). 
The superconductor has a magnetic penetration length A 
much larger than the coherence length £, and the film is 
taken to be a homogeneous thin plate of thickness d <C A. 
For the magnetic dots, we assume the magnetization on 
each dot to be quenched in, and pointing normal to the 
layer. In what follows, we will investigate how the inter- 
action between vortices in the superconductor is affected 
by the magnetic dots, and explore which quenched con- 
figurations of magnetization are favorable for vortex pin- 
ning at low temperatures. At the end, we will discuss 
how the desired configuration(s) of magnetization might 
be achieved in practice. 



FIG. 1. Schematic diagram of the superconducting layer 
(S) deposited over an array of magnetic dots (M), separated 
with a thin insulating layer (I). 

We first describe the case of a single dot with mag- 
netization M(r) to be specified shortly. Modeling the 
thin-film superconductor as an ideal sheet current K s (p) 
at the z = plane (with p = [x, y] denoting the position 
vector), Maxwell's equation of magnetostatics becomes 
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V 2 A(r) = —K s (p)5(z) 
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(1) 



in the gauge V • A = 0. We will describe the super- 
conductor in the London approximation which has 
K s (p) = (c/4ttA) ■ [5(p) - A(p )], where A = X 2 /d is the 
relevant magnetic length scale for the superconducting 
film, and $(p) is the London vector which, for a col- 
lection of vortices located at points p • with respective 
quantization (charges) tj, is 
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The above form of K s holds everywhere except at the vor- 
tex cores, which are normal regions of radius ~ £ around 
each singularity, and for as long as z ■ A(p ) = 0, which 
will be satisfied in our case. Using the superposition prin- 
ciple, we write A(r) = A s (r) + A m (r), where A s is the 
magnetic vector potential due to the supercurrent, and 
A m is the vector potential due to the magnetic dot. In 
what follows, we will model a magnetic dot by a perfect 
dipole of magnetic moment m = mi, placed at a distance 
£ below the superconducting plane and the origin. Thus, 
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A m (r) = m(z x r)/|r + £z\ 3 . The dipole model is ex- 
act, of course, only for homogeneous spherical magnetic 
dots. Nevertheless, we will use this dipolar approxima- 
tion throughout, in the hope that the results will give 
the correct order of magnitude for a variety of magnets 
whose magnetizations are normal to the plane. 

The free energy of the entire system may be written as 
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K s (p)| 2 + A s (p)-K s (p) 



— m • V x A, 



(3) 



The term in [. . .] is the sum of the kinetic and magnetic 
field energies of the supercurrent; the last term is the 
potential energy of the magnetic dipole in the magnetic 
field of the supercurrent. Using Eq. (|l|), we can integrate 
out the vector potential and express the free energy com- 
pletely in terms of the vortices in the superconductors. 
We have F = F vv + F vm , where 
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is the vortex- vortex interaction energy (including the vor- 
tex self-energy), with U{p) being the Pearl potential |t]] 
whose asymptotic behavior is 



U(p) 



ln(A/0, P « £, 

(l/2)ln(A/p), e«P«A, 
A/(2/>), p » A, 



and 
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E^(i^i), 



(5) 



(6) 



is the vortex-magnet interaction, with 

V(p) = - dn Mkp/A) 
Jo 2k + 1 

-A/21, p < £ <C A, 
-A/2p, £ < p < A, 
-2(A/p) 3 , p>A. 



(7) 



We will assume that p <C A and £ ~ £, which cor- 
respond to the more interesting and experimentally re- 
alistic cases. A single magnetic dipole with m > is in 
general able to bind more than one vortex. The attractive 
force on a vortex due to the magnet is the strongest at 
a distance ~ I from the magnet. Since the vortex- vortex 
repulsion decays only logarithmically with distance, we 
conclude that all bound vortices are concentrated in an 
area of radius ~ t around the dot. When their number 
is large (and it can be), they can be considered a single 
multiply quantized vortex, since £ ~ £ by assumption. 

The maximum charge n s of the bound vortex can be 
found as follows: Suppose the net vorticity of the sample 



(as enforced by the application of an external magnetic 
field) is AT > n s , and out of these N vortices, there is 
a vortex of charge n at the origin above the dot, with 
the remaining N — n vortices singly quantized and far 
removed from each other and from the origin. The total 
free energy of such a system is simply 



F(n) 
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(„2 + JV-n)ln^ + ^SnV(O). (8) 



n s can then be identified as the value of n which mini- 
mizes F(n), with the result 



u, = hit 
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where int[x] denotes the nearest integer to x. The numer- 
ical value of n s , which we shall call hereafter the satura- 
tion number, may be much larger than 1. As an example, 
let us consider a "typical" magnetic dot made of Co/Pt 
multilayer, with magnetization ~ 500 emu cm -3 0, dot 
size (0.25/xm) 2 and height h <~ 40 nm. Taking I to be 
h/2, we get m/£ ~ 3c/)q. The factor ln(A/£) is given by 
the superconducting film itself; for Nb thin films some- 
what below T c with £ ~ 20 nm, A/£ ~ 15 and thick- 
ness d ~ £ one has ln(A/£) ~ 5. In this case n s « 8. 
This example illustrates that magnetic pinning may be 
many times more effective than pinning by material de- 
fects (e.g., holes) at the same density, since the holes can 
directly bind only one vortex per site. The ability of a 
single magnet to bind so many vortices is perhaps the 
most desirable property of the magnetic dot system, at 
least for the purpose of vortex pinning; this property has 
however not been recognized previously [pj. 

Extrapolating the above-described properties of a sin- 
gle dot to systems with many dots, one might naively 
conclude that an array of dots in the ferromagnetic con- 
figuration would provide the strongest pinning. Indeed, if 
the applied magnetic field is n s times the matching field 
-B^, then each dot will bind n s field- induced vortices, a 
task not achievable by an array of non-magnetic pins if 
n s is large. However, the ferromagnetic dot array suf- 
fers from a different problem: It has difficulty pinning at 
small applied fields, i.e., for B <C n s B^, because of the 
appearance of intrinsic vortices and anti- vortices created 
by the magnetic field of the magnets themselves. 

To understand the effect associated with intrinsic vor- 
tices, let us again consider the single-dot system, but now 
in the absence of any external magnetic field. For weak 
magnets, no vortex is created in the superconducting 
film. Vortices eventually appear in the film for sufficiently 
strong magnets. This first happens when there is a doubly 
quantized vortex bound directly above the magnet, with 
two single anti- vortices straddling it, each at a distance po 
away Q . The anti- vortices must be present in a large film 
(of linear size ^> £), since the net magnetic flux through 
the x-y plane due to the magnet is zero. (This important 
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fact was left out of the model in Ref. f|.) Annihilation 
of the anti-vortices with the nucleus is prevented by the 
short-range repulsion between the anti-vortices and the 
magnet. To find the onset of intrinsic vortices, we com- 
pute the free energy F{p) obtained by applying the forms 
of the vortex- vortex and vortex-magnet interactions [Eqs. 
(|J) and (Q)] to the vortex-antivortex configuration de- 
scribed above. We take po to be the value of p which min- 
imizes F(p), with the result po = 167rm/(3(/>o). Intrinsic 
vortices appear when -F(po) < 0, since the free energy 
of the system without vortices is set at F — by defi- 
nition. In the experimentally relevant limit po 3> £j this 
occurs when the magnetization reaches a critical value 
m c = ^0o^ln(A/^). Using this expression for m c , we 
can rewrite Eq. QOj) as n s = int[| + ^-]- Note that 
n s = 2 when m — m c . Thus, intrinsic vortices and anti- 
vortices will appear once the magnet becomes a better 
pinning site than a simple void, i.e., for n s > 1. 

Next, let us consider the behavior of a ferromagnetic 
(FM) square array of dots with m > m c and lattice 
constant a ~ po, again in the absence of any external 
magnetic field. In the Co/Pt example mentioned above, 
a ~ po ~ 1pm. In this case, the intrinsic anti-vortices are 
very loosely associated with the magnetic dots and form 
a classical plasma, which interacts with the FM array 
only interstitially. For n s ^ 1, the interstitial pinning is 
of high order and therefore very weak. Consequently, the 
anti-vortices can be set into motion by a small applied 
current or thermal excitation, leading to dissipation even 
in the absence of any applied field! This is certainly not 
a desirable feature for the purpose of vortex pinning. 

The pinning ability of the FM array increases for in- 
creasing external fields, whose effect is mainly to annihi- 
late a fraction of the intrinsic anti- vortices at interstitial 
sites. The annihilation is complete when the applied field 
reaches the order of n s B$ as already described above. 
Thus maximum critical current is obtained at this field. 
Through this consideration, we see that the main effect 
produced by the ferromagnetic array is to shift the zero of 
the magnetic field to some large value, i.e., n s B^,. This 
can be utilized in special applications which requires a 
high critical current for a given field, but is not good in 
general, where high critical current is demanded for a 
range of external fields. 

The pinning properties of the sample in low applied 
fields can be significantly improved if the magnets form 
a quenched square antiferromagnetic (AFM) array, where 
the intrinsic anti-vortices produced by those magnets in 
the +z (or "up") direction tend to be attracted to the 
magnets in the — z (or "down") direction. For strong 
magnets whose binding distance pa (between the mag- 
net and its satellite anti-vortices) is of the order a, the 
anti-vortices will be pinned strongly to magnets pointed 
in the opposite direction, leading to a much larger criti- 
cal current needed to dislodge the intrinsic vortices and 
anti- vortices at low fields. 



In the absence of any external field, the free energy 
per unit cell of an AFM square array is obtained from 
Eqs. (||) and (||), assuming that n intrinsic vortices are 
localized on each upward magnets and n intrinsic anti- 
vortices are localized on each downward magnet. To sum 
the resulting alternating series of potential energies, a 
good approximation is to take the two largest terms of 
the series, which has a minimum when n = n s . 

In the presence of an external magnetic field, the field- 
induced vortices are subject to a periodic potential land- 
scape created by the magnets and the bound intrinsic 
vortices and anti-vortices. To characterize this potential, 
we place a "test" vortex of charge +1 at some position 
p, and compute the free energy W(p) experienced by 
this test vortex, due to interaction with all the intrinsic 
vortices and anti-vortices, riij = (— iy +3 n s , and with all 
the magnets, mij — (— l) J+J 'm, on the respective lattice 
sites Rij = a(ix + jy). The resulting expression 

is plotted in Fig. 2 for the example in the text with n s = 8 
and a = po/2. 




FIG. 2. Free energy W(p) of a test vortex in an AFM 
dot array. The range of W(p) is ±<jf>on s /(167r 2 A); p spans 
the primitive cell of the AFM lattice. The arrows indicate 
magnetic moments of the dots on the lattice sites. 

The form of W(p) clearly indicates that the test vor- 
tex will be attracted towards the nearest "up" magnet. 
Note that even though a magnet is already saturated 
with intrinsic vortices, it can still bind more field- induced 
vortices because there is no core energy cost associated 
with the latter. The maximum number of field-induced 
vortices n c that an upward magnet can bind is esti- 
mated as follows: Every field-induced vortex (per primi- 
tive cell) that is attracted to the upward magnet changes 
the potential landscape for the test vortex. These bound 
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vortices surround the nucleus, effectively increasing its 
charge. A test vortex will no longer tend to the magnet 
when the repulsion from the effective charge of the nu- 
cleus becomes equal to the maximum attractive force due 
to the magnet. The latter force is strongest at a finite 
distance D; in the point dipole approximation, Eq. (||), 
D ~ I. More realistically, D will be of order of the radius 
of the magnetic dot. The balance of forces gives 




2£ A 
— In. — — 1 
D £ 



(11) 



In our example system, where D ~ 0.1 /im, n c ~ n s . 

The analysis described in this study suggest very dif- 
ferent behaviors for a superconducting thin film on a 
(quenched) FM or AFM dot array. For the FM array, 
we expect the system to have low critical current J c at 
low fields, with successively increasing J c when the ap- 
plied field is an integer multiple of the matching field, 
and with the maximal J c obtained at B — n s B^. For the 
AFM array, we expect the pinning to be the strongest 
at low fields, with successively lower J c as the applied 
field reaches higher and higher orders of the matching 
field, and with the main effect diminishing beyond a field 
value of n c B,p/2. Thus, for different applications, one 
might want to use one or the other type of arrays, or 
some combination. 

It remains to address how the quenched magnetic dot 
configurations can be achieved and maintained in an ap- 
plied field. The FM array is most easily prepared by 
aligning the magnetic moments in a strong magnetic 
field in the z direction prior to measurement. It is less 
straightforward to ensure the AFM arrangement. Fortu- 
nately, we are aided in this case by the fact that the AFM 
state is the ground state of a square lattice of magnetic 
dipoles as long as the out-of-plane direction (i.e., the z 
axis) is the "easy" magnetization axis |flol] . The latter 
can be arranged by the construction of the individual 
magnetic dot, e.g., by using the multi-layer Co/Pt struc- 
ture |Q| . The appearance of intrinsic vortices further 
stabilizes the AFM structure. Thus, the AFM array will 
in fact form spontaneously upon cooling from high tem- 
perature in zero applied field, based on purely energetic 
considerations. There are, of course, kinetic constraints, 
such as coercive effects and the mobility of the domain 
wall separating the two degenerate states of the antifer- 
romagnet, that may prevent a perfect antiferromagnet 
from being formed at a reasonable time. These kinetic 
constraints are, on the other hand, necessary to prevent 
the magnets to undergo a spin-flop type transition |T^ ] 
to the FM phase upon increasing the external magnetic 
field. One can also envision building an array of micro- 
scopic superconducting current rings below the film, and 
electronically control the sense of current in each ring. 
This way, arbitrary quenched magnet configurations can 
be specified. 



FIG. 3. Schematic diagram of the superconducting layer 
with thin magnetic rods magnetized horizontally. 

A practical way to achieve effects similar to that of the 
AFM array is to bury under the superconductor an ar- 
ray of thin magnetic bars (Fig. 3). The shape anisotropy 
of the bar will force the magnetization to be along the 
length of the magnet. The field outside the magnet is 
the same as that of two oppositely charged magnetic 
monopoles located at its endpoints, creating an effect 
similar to that of an antiferromagnetic pair of dipoles 
whose moments are normal to the sample p3[ . A crude 
approximation to this geometry is the use of thin mag- 
netic disks whose diameter is a significant fraction of the 
lattice constant. The anisotropic shape of the disk in- 
duces a magnetization parallel to the superconducting 
plane, similar to that of the magnetic bar. This was in 
fact the geometry used by Martin et al. [0. 
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